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We investigate the Yukawa model in which Nf fermions are coupled with a scalar field <j) through 
a Yukawa interaction. The phase diagram is rather well understood. If the fermions are massless, 
there is a chiral transition at T c : for T < T c chiral symmetry is spontaneously broken. At Nf = 
°° the transition is mean-field like, while, for any finite Nf, standard arguments predict Ising 
behavior. This apparent contradiction has been explained by Kogut et al., who showed by scaling 
arguments and Monte Carlo simulations that in the large-rV/ limit the width of the Ising critical 
region scales as a power of 1 /Nf, so that only mean-field behavior is observed for Nf strictly 
equal to infinity. We will show how the results of Kogut et al. can be recovered analytically in the 
framework of a generalized l/Nf expansion. The method we use is a simple generalization of the 
method we have recently applied to a two-dimensional generalized Heisenberg model. 
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Finite-temperature transitions in quantum field theory models have been the object of many 
theoretical studies. Here we investigate the transition in the 2+1 Yukawa model (but the arguments 
can be generalized to generic (d+l) models, d < 4) in which a scalar field is coupled to Nf de- 
generate fermions by a Yukawa interaction. As discussed in Ref. [jl|] this model shows a peculiar 
behavior for Nf — > °°. For finite Nf, dimensional reduction predicts that the finite-temperature tran- 
sition, if continuous, belongs to the two-dimensional Ising universality class. On the other hand, 
for Nf = °° an explicit calculation gives mean-field behavior [0]. These two apparently contradic- 
tory results were explained in ^ in terms of a critical-region suppression. A similar behavior was 
observed in a generalized 0(N) G model in [J3j] and an explanation was provided in [Q]. The same 
techniques developed in [Q] can be applied here to obtain an analytic description of the crossover 
from mean-field to Ising behavior in the large-iV/ limit. 

We will investigate the Yukawa model with action ^ 



I" JVy 

S>[4>,Y,Y} = Jd 3 r ^(d<i>) 2 +Nf^<i> 2 +Nf-<i> 4 + pWf(d+g<!>+M)y f 



(1) 



where the integral is over M 2 x [0, 1/7] and we use periodic (antiperiodic) boundary conditions for 
the boson (fermion) in the "temporal" direction. We imagine the theory to be somehow regularized 
with a cutoff that sets the energy scale. The chosen regularization is not relevant for the discussion, 
and, for instance, the reader may imagine using the lattice action with staggered fermions studied 
in©. 

In order to determine the behavior for Nf = °° we integrate out the fermionic degrees of free- 
dom. Starting from ([I]) we obtain 

r N >' 

e -N f y eS m = / Y\ dl ff fd¥fe -yi^vM (2) 

&tsm = ^{d^) 2 + ^ 2 + -^-Xr\og{i+g^+M) . (3) 

For Nf = oo we can perform an expansion around a translation-invariant saddle point . The sta- 
tionarity condition gives the gap equation 

— f A d 2 p 1 

am + ^Xm 3 = (m + M)T V / . . _ -= — — ^ — (4) 

6 f&J (2n) 2 p 2 + (m + M) 2 + (2n+l) 2 n 2 T 2 



Here we have defined m = g<p,~jl = p.g~ 2 , and A = Xg~ 4 . With this choice, the dependence on g 
disappears. Note that the gap equation is symmetric under m — > —m, and M — ► — M so that it is not 
restrictive to consider m > 0. 

We first analyze the model for T = 0. A simple analysis of the gap equation shows that there 
are two possibilities for M = 0. There is a critical value ~jl c such that for ~jl <Jl c the relevant 
solution of the gap equation is such that m^O (chiral symmetry is broken), while in the opposite 
case we have m = (no chiral symmetry breaking). For the behavior is smooth in the 

parameters. In the following we will only be interested in the case in which the zero-temperature 
theory shows chiral symmetry breaking. Thus, we shall assume JX < JI C . Parameters A and JI do 
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not play any additional role in the model and thus in the following we will not consider explicitly 
the dependence on them. 

Let us now consider the behavior for T ^ 0. We restrict the discussion to ~jX > 0; in this case 
symmetry is restored for T — > °° and M = and therefore, there exists a value T c such that for 
T < T c , m is a nonvanishing function of T (chiral symmetry is broken), while for T > T c , m = 
(chiral symmetry is restored). Explicitly, by using (Q) we obtain the relation 

-_ T y [ A jfj> 1 

If M 0, the behavior is smooth in T. Thus, T = T c , M = is a critical point (CP) for model ([]]). 

Close to the CP one can define a thermal scaling field u t and a magnetic scaling field uj,. Using 
the gap equation it is easy to see that we can take u t = (T — T c )/T c and «/, = M/T c . Then, for 

—2/3 

u t ,ut, — > at fixed x = u t u h , one obtains the mean-field equation of state 

m/T c = u h 1/3 f(x) af(x) 3 + bxf(x) + 1=0 (6) 

with a > 0. From (^), if u t = 0, we obtain m ~ w/, 1 / 3 ', so that 8 = 5mf = 3; if Uh = 0, we have m ~ 
M; 1 / 2 so that j6 = j8mf = 1 /2. Thus, for Afy- = °°, the behavior is of mean-field type, in agreement 
with previous results [0]. Note that the condensate £ = (i/fi//) is proportional to M + m « ui^ 3 f(x) 
in the scaling limit. A completely analogous discussion holds for the staggered lattice model. 

In order to perform the l/Nf calculation, we expand action <j3j> around the saddle-point solu- 
tion, writing = + 0/ <JWf . We obtain the expansion 



= \Lj 70 i 0, 1 (p)a- 1 (p)0-«(-p) 



0n(p) 0m (q) 0-„- w (-P ~ q) 

1 /■ d 2 p d 2 q d 2 k - f4 x , 

+ 4!a^,/ ^;q,-;M;-P-q-k,-m-n-,) 



(2tt) 2 (2tt) 2 (2tt) 2 

0n(p) 0m (q) 0f (k) 0-«-m-r(-p - q - k), (7) 



3 /2 

where the neglected terms are of order l/NJ . Here 0„(p) is the Fourier transform of 0(*,r) 

(x £ [0,r'],r,peR 2 ): 

0(x,r) = £ e 2 ™ T f ;1 (p)e'P r . (8) 

neZ J \ zn > 

In the standard approach, 1 /Nf expansions are obtained by performing a perturbative expansion of 
theory (0) in powers of 1 /Nf. This is possible here only far from the CP, since at the CP the field 0o 
is massless and the expansion is plagued by infrared divergences. Indeed, starting from the explicit 
expression 

-2a-1/ an -, l T2 f A d 2 q q 2 + (2n+l) 2 K 2 T 2 -(m + M) 2 

g A n (P = 0) = LL + -Am Z - T > / 7^5 ^ ,„ 7 77^777 (9) 

° y 1 ^2 £ Z J (2n) 2 [q 2 + {2n+l) 2 K 2 T 2 + {m+M) 2 } 2 
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and using (Q), we find that close to the CP 

A l (0) ~ m 2 ,(m + M) 2 (10) 

No singularity arises for the other modes 0„, since A~ x (p = 0) = (iKtiT) 2 at the CP. The infrared 
singularity of Ao(p) is of course expected: A,7 ! (p) is proportional to the square of the mass of the 
boson field that should vanish at the CP (at the CP the correlation length diverges). 

In order to deal with this singularity we use the technique we have recently applied in [Q] to 
a generalized Heisenberg model in two dimensions. We first integrate out the massive modes and 
obtain an effective action for the zero mode: 

e -y,^\ = [Hdbe-^M. (11) 
The effective action ^ e ff[<p] has an expansion in powers of 1/Nf of the form 



y eS [(p] = //<p(0) + iyV(p)+r]<p(p)<p(-p 



p 



3! 



/ / v J (p ) q>-p-q) < p(p) < p(q) < p(-p-q) 



pjq 



+^J^JJv {4) (p,q,s,-p-q-s)(p(p)q)(q)(p{s)(p{-p-q-s), (12) 

where <jp = a<j>o + b,u = c/Nf, and a, b, c are functions of T and M fixed by the following normal- 
ization conditions: 

K(p) = p 2 + <9(p 4 ) for p -» 

(0,0, 0,0) = 1 (13) 
7^(0,0,0) =0 

for any T and M close to the CP. The last condition is not trivial and can be imposed because the 
same property holds at the CP for V® (0,0; 0,0; 0,0). A gain, in ( |12| ) we have neglected higher- 
order terms in 1/Nf. 

At this point, the origin of the mean-field-to-Ising crossover is quite clear. For Nf = °°(u = 0) 
the zero mode is a free field and thus the model shows mean-field behavior. On the other hand, 
for finite Nf, one must consider the full theory (ff2[), which is nothing but a generalized <p 4 theory 
whose critical behavior belongs to the Ising universality class. Note that in (|l2|) we have discarded 
higher-order vertices <p" that are multiplied by higher powers of 1/Nf: one can show that they do 
not play any role in the crossover limit we shall describe below ffi\. 

We wish now to compute the crossover behavior. Formally, the momentum dependence of the 
vertices is irrelevant (one may think of the formal continuum limit in lattice theories) and thus we 
can simply consider (this relation becomes exact in the infrared limit, in which one only considers 
the long-distance behavior; for a proof in a lattice framework, see 

^cont[<p] = / d 2 r [H(p + \(d(p) 2 + i_r(p 2 + ±,u(p 4 ] . (14) 

Then, we define %(r) = <p(r/ y/u) and note that we can rewrite the action as 



^contM = J d 2 r 



H 1 r ? 1a 

u 2 2u 4! 



(15) 
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In terms of x, the action is a function of H/u and r/u. Then, consider the zero-momentum con- 
nected correlation function % n . We have 

Xn = J d 2 r 2 . ..d 2 r n (<p(0)<p(r 2 ) . . . «p(r„)) conn 

J d 2 r 2 ...d 2 r n (x(0)x(r 2 )...x(rn)) conn = u n - 1 MH/u,r/u), (16) 



i.e. u l ~"Xn is a scaling function of H/u and r/u. Unfortunately, the derivation is not correct, since 
we have not taken into account the presence of the cutoff that breaks scale invariance. However, 
in two dimensions (and, in general, for d < 4) only a mass renormalization (a redefinition of the 
parameter r) is needed in order to take care of divergencies. By a proper treatment [||] one can 
show that there is a function r c (u) (in two dimensions the determination of this function requires 
only a one-loop computation since the only diverging diagram is the tadpole) such that, for t = 
r — r c (u) — ► (infrared limit), H — ► 0, and u — > (weak-coupling limit), the correlation function Xn 



satisfies the scaling relation Xn = u"~ l f„(H/u,t/u). These results extend to theory (12) although 
the presence of odd powers of (p requires also a renormalization of H. In we showed that one 
can find functions r c (u) and H c (u), such that, for t = r — r c (u) — > 0, h = H — H c (u) — > 0, u — > at 
fixed h/u and t/u, the correlation function Xn scales as 

Xn~U ] -"f n (h/u,t/u). (17) 

Functions f n (x,y) are universal: they do not depend on the explict form of the vertices and of K(p) 
and can be computed directly in the continuum theory. They are the crossover functions that relate 
mean-field and Ising behavior. Consider, for instance, the case h = 0. For t fixed and i/-»0we 
obtain the standard perturbative expansion; thus, t/u — > oo corresponds to the mean-field limit. On 
the other hand, for t — > at u fixed, Ising behavior is obtained; t/u = is the nonclassical limit. By 
varying t/u between and oo one obtains the full universal crossover behavior. The universality of 
the crossover implies that these functions can be computed in completely different settings: one can 
use field theory [||, ||], generalized Heisenberg models such as those considered in [Q], or medium- 
range models (see [||, [7]] and references therein). For instance, crossover curves for the effective 
exponents for M = can be found in [Q] (field theory) and in ^ (in this case the correspondence 
is R 2 = bNf, where the constant b can be computed by matching the corresponding perturbative 
expansions at one loop). 

The analysis of [Q] can be simplified in the present model since chiral symmetry is preserved 
by regularization. In this case the relations 

H = V (3) (p,q,r) =0 (18) 

hold at the CP. They imply H c (u) = 0. 

In theory ( ff2[ ) r, H, r c (u), and u are functions of T, M, Nf. Thus, the next step consists in 
determining how these quantities should scale close to the CP in order to keep x t =Nf{r — r c (Nf)) 
and x h = Nf(H — H c (Nf)) fixed. A calculation gives [Q] 

a + b\oe,Nf 

T c (N f ) - T c (Nf = oo) « ^JL^, (19) 
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where a and b can be explicitly computed. Moreover, we have 

T-T c (N f ) x t M 



x h 



T e (N f ) 



T c (N f ) Nf^ 2 ' 



(20) 



Ising behavior is observed for x t ,Xh <C 1. This confirms the critical-region suppression predicted 
in ^ : The width of the Ising critical region scales as NJ 1 in the thermal direction and as 3/<2 in 
the magnetic one. For M = one can also characterize the crossover in terms of the bosonic mass 
Mbos- Indeed, similar arguments (see [FJ]) allow us to predict 



N 



m f M [x t =N f (T-T c )/T c 



(21) 



/ 



where /m(x) is a crossover function that behaves as x for x 1 (V = 1 in the Ising theory) and 
as *Jx for x 3> 1 (v = 1/2 for a Gaussian field). This relation shows that fixing x t is essentially 

1/2 

equivalent to fixing xm = Ay M\, os /T c : Ising behavior is observed for xm -C 1, mean-field behavior 
in the opposite case. This is exactly the scaling condition discussed in [jlj]. 
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